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ABSTRACT 

In the symbiotic (or WD+RG) channel of the single degenerate scenario for 
type la supernovae (SNe la) the explosions occur a relatively long time after star 
formation. The birthrate from this channel would be too low to account for all 
observed SNe la were it not for some mechanism to enhance the rate of accretion 
on to the white dwarf. A tidally enhanced stellar wind, of the type which has 
been postulated to explain many phenomena related to giant star evolution in 
binary systems, can do this. Compared to mass-stripping this model extends 
the space of SNe la progenitors to longer orbital periods and hence increases 
the birthrate to about 0.0069 yr -1 for the symbiotic channel. Two symbiotic 
stars, T CrB and RS Oph, considered to be the most likely progenitors of SNe la 
through the symbiotic channel, are well inside the period-companion mass space 
predicted by our models. 

Subject headings: binaries: close — stars: evolution — supernovae: general 



Introduction 



Type la supernovae (SNe la) are used as distance indicators for cosmology. From 



SNe la we deduce tha t the expansion of the Universe is accelerating (IRiess et al.l Il998 



Perlmutter et al.lll999l ). They are also crucial for the study of Galactic chemical evolution 
because they expel iron. However, the exact nature of their progenitors remains unclear and 
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it is still uncertain whether a double-degenerate (DD) or single-degenerate (SD) scenario 
dominates. In the DD scenario two carbon-oxygen (CO) white dwarfs (WDs) can produce 
a SN la while merging if their total mass is larger than 1.38 M & . In the SD scenario, a CO 
WD explodes as a SN la if its mass reaches 1.38 M & while accreting from a non-degenerate 
companion such as a main-sequence (MS) star or slightly evolved subgiant in the WD+MS 
channel, a giant star in the WD+RG channel or the s ymbiotic channel or a helium star 
(Ivan den Heuvel et al.lll992l ; Irian k Podsiadlowskil 120041 ) . 



A progenitor model can be tested by comparing the predicted distribution of the time 
between star formation and SN la explosion with that observed. This ranges from less than 
10 8 yr to more than 3 Gyr. A significant number of SNe la have taken the longest times 
( jBotticella et al.ll2008b iTotani et alll2008uSchawinskill2009l ). Thou gh the DD scenario can re- 
produ ce such delays the SD scenario cannot, except in the study by lHachisu. Kato &: Nomoto 
( 120081 ) who assumed that the wind from the WD strips the outer layers of the red giant (RG) 
at a high rate (their mass-stripping effect). They showed that s ymbiotic stars are likely pro- 
genitors of SNe la with long delays. iGilfanov fc Bogdanl (120101 ) have claimed that supersoft 
X-ray (SSX) fluxes in early type galaxies are m uch smaller than expected i f the SD scenario 
were to dominate SNe la production. However lHachisu. Kato fc Nomotol ( 120101 ) argue that 
this SSX flux is actually rather strong evidence in favour of the SD scenario in early type 
galaxies and that the progenitors are symbiotic stars. 

There are at least four recurrent symbiotic novae, such as RS Oph or T CrB, that have 
a very massive WDs. They are t hought to have red giant companions of mass M < 1 M & 
( lAnupama fc Mikolajewskal 119991 ). It is difficult to reproduce these systems in standard 
binary evolu tion because dynamical ins t ability in mass transfer sets in when the giant fills its 
Roche lobe. lHachisu. Kato fc Nomotol ( 1l999a[ ) systematically studied binary evolution and 
established WD+MS and WD+RG paths to SNe la by introducing three new physical effects, 
(i) the wind from the RG acts like a common envelope to reduce the separation for very wide 
binaries with separations up to about 30,000 R Q , (ii) the WD loses much of the transferred 
mass in a massive optically thick wind and (iii) their mass-stripping effect. They produced 
SNe la at a rate comparable with that observed. Since then the only work to produce 
sufficiently high SN la rates with long delay times was Hachisu, Kato &; Nomoto 
(2008), who implemented all the three effects. Other works adopted the above 
effects except mass-stripping in their population syntheses and found smaller SN 
la rates. 

The concept of a compact object accreting from the stellar wind of a companion was used 
by iDavidson fe Ostrikerl ( 119731 ) to explain X-ray pulses in massive X-r ay binaries. It is now 
generally included in population synthesis calculations such as those by lHurley. Tout fc Pols 
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( 120021 ). though lHachisu. Kato &: Nomotol (ll999al ) did not include any wind accretion in their 
model. The surfaces of red giant stars are not tightly bound and so it is relatively easy to drive 
a cool wind. Such winds might be enhanced by tidal or other interaction with a companion 
( ITout & Eggletonlll988al ) and so could remove significant mass an d angular momentum from 
the giant's envelope before Roche lobe overflow (RLOF) begins. ITout fc Eggletonl ( 1l988af ) 
introduced this concept of companion reinforced attrition (CRAP) to explain mass inversion 
in RS CVn binaries and it has since been widely used to explain phenomena related to 
giant star evolution in binary systems, such as double degenerate objects, post-AGB s tars, 
barium stars, cataclysmic variable stars and bipolar plan etary nebulae (see for example lHan 



19981 ; Ivan Winckelll2 003; Bo nacic. G lcbbce k &: Polsll2008l ) . As potential progenitors of SNe la 



symbiotic stars ( IKenyon &: Webbink! 119841 ) gain two advantages from CRAP. First the WD 
may grow in mass substantially by accretion from stellar wind before RLOF and secondly 
mass transfer may be stabilized because the mass ratio (M giant /M WD ) can be much reduced 
at the onset of RLOF. So the distribution of masses and periods from which SNe la arise 
through the symbiotic channel is enlarg ed when CRAP is taken int o accou nt. Here we study 
this for SNe la with CRAP instead of lHachisu. Kato fc Nomotol (Il999af )'s mass stripping 
effect. Although both processes may occur in real systems we include only CRAP so as to 
isolate its effects. 



2. Binary Evolution Calculations and Results 



We evolve about 2,000 population I (wit h metallicity Z = 0.02) close WD +MS binary 
stars with Eggleton's stellar evolution code ( lEggletonl Il971t iPols et al.l 119981 ). Their CO 
WD masses Mwd axe 0.6, 0.65, 0.8, 1.0 and 1.1 M Q , the initial orbital periods are between 
10 and 10,000 d spaced at intervals of Alog 10 (Pyd) = 0.1 and the companion, secondary 
masses M 2 in the range from 0.8 to 7.0 M G) . T he tidally enhanced mass- loss rate from the 



sec ondary Mi^ is m o delled by iReimersi (119751 )' formula with the extra tidal term included 
bv lTout fc Eggletonl Jl988bh so that 



M, 



2w 
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-13 



0; 
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1 R_ 
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(1) 



where R and L are the radius and luminosity of the giant secondary, is its Roche lobe 
radius, 77 is Refiners' wind coefficient. We set 77 = 0.25. The wind enh ancement parameter 
B w we set to 10,000 as required to fit Z Her (ITout fc Eggletonl Il988al ). Thus the mass- loss 
rate |M 2w | could be 150 times as larger than Reimers' rate when the star more than half fills 
its Roche lobe. 
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Some of the mass lost in the stellar wind of the giant may be accreted by the WD at a 
rate (jBoffin &: Jorissenlll988l ) so that 



M' 



1 



2a 



vT 



_ ( GM WD 



2\v 



(2) 



where t> or b = a/ G(M 2 + Mwd) / a is the orbital velocity, G is Newton's gravitational constant, 
a is the semi-major axis of the orbit and e is its eccentricity We take e = because we 
expect orbits to have circularized by this time. The coefficient ct acc is an accretion efficiency 
and we set a acc = 1.5. We fix the wind velocity v w to Skms" 1 and discuss this later. If 
v w or a is small, the right hand side of equation (j2J) becomes large and so we add the limit 
A^2a < —M 2w . That part of the wind which is not accreted carries off the specific angular 
momentum of the donor star. 

Before RLOF begins stellar wind is the only way to transfer material to the WD so the 
mass transfer rate, M t = M 2a . During RLOF the mass is transferred in both a stream and 
the wind so that M t = M 2a + |M 2 rlof|, where M 2 rlof is mass-transfer rate by RLOF. 



We continue the practice of limiting the rate a t which the WD can a ccrete (|Han fc Podsiadlowski 



2004! ; iMeng. Chen fc Hanll2009l ; IWang et al.ll2009f ) by the prescription of lHachisu et all fll999bh 
This accounts for the limited rate at which hydrogen can burn and for mass loss in both 
hydrogen and helium novae. Thus 



M WD = r] R r] He M t . 



(3) 



Hydrogen accretion is controlled by 



M C1 
| Aft 

1, 

o. 



\M t \> M cr , 

Mcr > \M t \ > M low , 
\M t \ < M low , 



(4) 



where M\ ow , equal to \M cn is the accretion rate below which hydrogen novae expel most of 
the material and M cr is the critical accretion rate above which hydrogen cannot burn as it is 
accreted but is instead expelled in an optically thick wind. Helium accretion is then further 
controlled by 



-O.^ilog^M^/MQji- 1 ) + 5.35) 2 + 1.05, -7.3 < log 10 (M Hc /M yr" 1 ) < -5.9, 

-5.9 < log lo (M He /M yr- 1 ) < -5, 

where Mn e = ?7h I M t \ . That part of the mass transferred but ultimately not accreted carries 
off the specific orbital angular momentum of the white dwarf. 
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In Figure [T] we show the evolution of a binary system which ends as a SN la. Ini- 
tially = 0.8 M , M\ = 1.5 M and log 10 (P7 d ) = !- 9 - Tne wind enhancement factor 
£>w = 10,000. We see three phases of CO WD growth to 1.378 M at which carbon ig- 
nites degenerately. First, wind accretion becomes important at about 4.7 x 10 6 yr when the 
accretion is sufficient for hydrogen to burn as it accretes. Stable mass transfer dominates 
for a short time of around 1.02 x 10 7 yr. Finally wind accretion continues after the system 
has detached because the orbit grows faster than the giant. The supernova occurs when the 
secondary's mass has fallen to 0.353 M just as it is about to leave the red giant branch and 
shrink to a white dwarf. 

At 2.673 x 10 9 yr (t ~ 7 x 10 6 yr in Figured]) the enhanced mass-loss rate reaches its 
maximum of 150 times Reimers' rate. In this case the WD is able to accrete all the wind 
and M 2a = —M 2w . This is an extreme case of small P 1 , corresponding to the short orbital 
period edge of the space that leads to SNe la for an initially 0.8 M & CO WD (see Figure [2] 
below). At longer initial periods the WD reaches 1.378 M before RLOF occurs so there is 
no stable mass transfer phase. For much wider systems and for more massive secondaries, 
that ignite helium non-degenerately and consequently do not grow very large on their first 
ascent of the giant branch, similar evolution occurs when the secondary is on the asymptotic 
giant branch. 

Figure [2] shows the results of our binary evolution calculations in an initial period- 
secondary mass plane for various initial WD masses. No SNe la arise when = 0.6 M . 
Whether a CO WD reaches the critical mass to explode as a SN la depends on both the 
amount of mass transferred and the efficiency at which this can be accreted by the CO WD. 
In each panel contours enclose the progenitors of SNe la. Only for the very high initial 
white dwarf masses, of 1.1 M , can a massive secondary star drive a SN la. These stars 
ignite helium non-degenerately and so do not grow large enough to transfer sufficient mass 
while on their first ascent of the giant branch. Thus, only those initially massive enough to 
retain enough envelope mass to transfer via a wind on the first ascent of the giant branch 
are able to reach SNe la. In this case the initial mass of the secondary ought not to be larger 
than the progenitor of the WD. This determines our upper limit for the upper region in 
panel (d) and explains why this region is absent for lower-mass WDs. The other regions are 
more interesting because they are more populated by virtue of the stellar mass function. For 
these the low secondary mass boundary is due to the limited amount of mass available for 
transfer from the lower-mass companions. The long period boundary is due to the efficiency 
of burning and accumulating the transferred material. The high secondary mass and low 
period boundary are both due to the onset of RLOF and dynamically unstable mass transfer 
because CRAP has not removed enough mass. Carbon burning then begins more gently in 
the heated outer layers of the WD and a supernova is avoided. 
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To examine the effect of CRAP in more detail we repeated our calculations for an 
initially 1 M WD with < B w < 10,000 (see Figure [3]). As B w decreas es the progenitor 
region shrinks and vanishes as B w — > 0. This is consistent with the claim of llben &: Tutukov 
(119841 ) th at there is no symbiotic channel t o SNe la. For comparison, the progenitor region 
found by lHachisu. Kato fc Nomotd (11999al ) is also shown in this figure. Our model extends 
the possible progenitor periods by almost a factor of ~ 4 — 5. 

The velocity v w of the companion's wind is important because it strongly affects M 2a . 
However the growth rate of the CO WD M WD is not always sensitive to v w because M WD is 
limited by M cr when the mass accretion rate is large enough (equation H]). Thus once M 2w is 
larger than about 10 -6 M Q yr _1 the dependence of Mwd on v w becomes small because only 
about 10 percent or less of M 2w is required for the CO WD mass to increase at M cr . When 
M 2w is closer to 10~ 7 M yr _1 the CO WD may not increase at all if v w is large. Our choice 
of v w = Skms" 1 is an estimate of the lower limit to the poorly modelled wind velocity. We 
choose it to demonstrate that such a channel to SNe la can exist but note that our progenitor 
region would shrink in places if v w were largeiQ. 

Figure H] shows the regions in the final period-secondary mass space at which the WD 
reaches 1.378 M for various initial C O WD masses. The region for a 1 M CO WD cal- 
culated with the mass-stripping effect (lHachisu. Kato fc Nomotd 120081 ) is superimposed in 
the figure. As for the initial systems, these regions extend to longer orbital periods. Two 
symbiotic stars, that are considered possible progenitors of SNe la, T CrB and RS Oph, with 
orbital periods 228 d (IBelczyriski fc Mikolajewskalll998l ) and 454 d (IBrandi et all 120091 ). are 
plotted under the assumption that Mwd ~ 1-378 M . Both are located well inside the region 
of SNe la progenitor regions for our models but lie on the edge of what might be expected 
with mass stripping. 



3. Frequency estimates 

Figure [2] shows how the CRAP enhanced symbiotic channel would be very effective 
if many binary stars with a massive CO WD and an unevolved companion could form 
with periods between about 100 and 1,000 d. However conventional population synthesis 



1 We investigated two alternatives in a little more detail (a) v w = lOkms -1 and (b) v w — a w Woo, where 
Vnc is a terminal wind velocity of typically 15kms _1 for a giant star and a w is calculated according to 
Yungelson et a In both cases the progenitor regions shrink with the longest period systems not 



reaching sufficient WD masses to ignite. For example the upper period boundary moves to log 10 (P/d) > 2.9 
for M WD = 1.0 Mq and no SNe la are predicted at all when M^ VD = 0.65 M Q in case (a). 
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calculations tend to produce such systems either with P < 100 d following the ejection of 
a common envelope or with P > 1,000 d if they have avoided RLOF. On the other hand 
we do observe symbiot ic stars, such as T CrB and RS O ph, with massive WDs and orbital 
periods around lyr. lHachisu. Kato fc Nomotol (ll999at ) describe one way in which these 



systems, and consequently many like them, could have formed by shrinking a wide orbit in 
a process similar to common envelope evolution but by interaction with the stellar wind. 
To obtain the frequency of SNe la produced by our channel we really ought to carry out a 
full populatio n synthesis from the zero- a ge mai n sequence but here we apply the procedure 
explained by lHachisu. Kato &: Nomotol (Il999al ) in their sections 4.3 and 4.4 so as to get a 
direct comparison with their rate. We divide the WD masses into t h e ran ges given in our 
Table [T]and calculate the integral in equation (1) of llben fc Tutukovl (119841 ). 



0.2yr _1 AgA log 



10 ' 



Mr, 



M A 



M~ 2 - 5 dM, 



(6) 



where Ma and Mb are the typical progenitor masses of the WDs at each end of the range. 
The integral then accounts for the initial mass function. We calculate the initial range of 
mass ratios, Aq, for a given WD mass by 



Aq 



M u M 



M A 



M b ' 



(7) 



where M u and M\ are the upper and lower limits to the SNe la progenitor regions identified 
directly from Figure EJ By formula (|6]) we assume that the mass ratio is uniformly distributed 
between and 1. With the quantity A log 10 a we assume that the initial period distribution is 
flat in log P and that this difference is unchanged by the evolution up to the point of Figured] 
so that we may read it from the figure directly for the particular WD mass. We then calculate 
a simple mean rate v of the tw o extremes. Our total SNe la rate is 0.0069 yr _1 , significantly 
larger than that calculated by lHachisu. Kato fc Nomotol (Il999al ) with their mass-stripping. 
This rate is very uncertain and probably overestimated because we, as did they, have used 
cuboid shaped regions in the initial parameter space but it is not dissimilar to the actual 
observed rate for a galaxy like our own. 
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Table V. The frequency of SNe la from our model with B w = 
10,000. For each range of WD masses M WD , Ma and M B are the 
corresponding lower and upper zero-age progenitor masses and Aq 
and A log 10 a are the typical mass ratio and separation ranges for 
the regions in Figure [2] which lead to SNe la. The frequencies v 
and v are the calculated supernova rates if the extremes of the 
distribution apply across the whole range and a simple mean of 
these two extremes. The factors of 2/3 account for the change from 
period to separation space. 





A log 10 a 


M A /M 


M B /M Q 


Aq 


zz/yr" 1 


v/yr 1 


0.65-0.8 


0.3 x 


2 
3 


2.42 


4.48 


0.3168 


0.0013 


0.0043 




0.8 x 


2 
3 


2.42 


4.48 


0.6412 


0.0073 




0.8-1.0 


0.8 x 


2 
3 


4.48 


6.63 


0.3101 


0.0001 


0.0019 




1.7 x 


2 
3 


4.48 


6.63 


0.3960 


0.0028 




1.0-1.1 


1.7 x 


2 
3 


6.63 


7.58 


0.2376 


0.0004 


0.00055 


(lower contour) 


2.0 x 


2 
3 


6.63 


7.58 


0.3787 


0.0007 




1.0-1.1 













0.0000 


0.00015 


(upper contour) 


0.7 x 


2 
3 


6.63 


7.58 


0.4489 


0.0003 
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Fig. 1. — The evolution of a binary system that leads to a SN la at the point marked with 
the filled circle on the vertical line just after 1.2 x 10 7 yr, under the companion reinforced 
attrition process with B w = 10,000. The time is arbitrarily offset by 2.666 x 10 9 yr from the 
birth of the system so that the timescale for the accretion process is more easily seen in this 
figure. The evolution of the secondary mass M2, the CO WD mass Mwd, the mass loss rate 
from the secondary — M 2 , the mass-growth rate of the CO WD Mwd and the orbital period, 
log 10 (P/d) — 1.2 are shown. For this binary, which is an extreme case (see discussions in the 
text), the mass-transfer rate M t = |M 2 |. The peak at t « 1.02 x 10 7 yr marks a phase of 
stable RLOF. The CO WD growth in other models is similar but most have wind accretion 
only before RLOF. 
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Fig. 2. — Results of binary evolution calculations in (logP 1 , M^) planes, where B w = 10000 
for four initial masses M^ D . Initial periods are P 1 and initial secondary masses are M\. Filled 
symbols indicate the progenitors of SNe la explosions. For filled triangles and squares the 
explosion occurs before helium ignition in the secondary star. For the triangles the explosion 
occurs during RLOF while for the squares the explosion occurs in the wind accretion phase. 
The final companion remnant would be a helium white dwarf in these cases. Filled circles 
indicate systems that explode as SNe la after central helium burning and before RLOF. In 
these cases the final companion remnants would be CO WDs. Symbols ' x ' and stars indicate 
systems that are unstable to dynamical mass transfer with CO WD masses lower or higher 
than 1.35M respectively but still less than 1.378 M . Open circles are systems stable to 
dynamical mass transfer but in which there is insufficient mass for the CO WDs to reach 
1.378M . Systems marked '+' and ! ®' begin AGB thermal pulses before RLOF or a SN la. 
It becomes computationally difficult and inefficient to evolve these stars fully so we make 
a simple estimate of whether a SN la is likely. Only in the system marked is the CO 
WD likely to reach 1.378M . The heavy contours enclose the systems in which the WD can 
explode as a SN la. 
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Fig. 3. — Contours for SNe la in the (logP 1 , M^) plane for a 1 M & CO WD with various values 
of B w . The thin solid triangle roughly enclos es the SNe la progenitors in the mass-stripping 
scenario of lHachisu. Kato fc Nomotol (Il999al ). 





log 10 (P SN /d) 

Fig. 4. — Contours enclosing progenitors of SNe la, now in the final period-secondary 
mass plane just before explosion for various initial WD masses and also the triangle for 
Hachisu. Kato fc Nomotol (jl999af Ts mass-stripping scenario. Filled stars and circles are the 
current positions of two symbiotic binary stars, T CrB and RS Oph. To plot these we have 
assumed both the CO W Ds have a mass of 1.378 M m so tha t the secondary mass for T CrB 
is 0.71 M Q according to Hachisu. Kato fc Nomotol f 1999a) or 0.823 M m derived from the 
mass ratio of about 0.6 de termined by (IBelczyriski fc Mikolajewskalll998l ) and for RS Oph 
0.68 — O.8M according to iBrandi et al.l ( 120091 ). As in panel (d) of Figure [2]the contour for 
1.1 M & has two parts but the upper region is off this figure. 



